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We study magnon modes in the presence of a vortex in a circular easy-plane ferromagnet. The 
problem of vortex-magnon scattering is investigated for partial modes with different values of the 
azimuthal quantum number m over a wide range of wave numbers. The analysis was done by com- 
bining analytical and numerical calculations in the continuum limit with numerical diagonalization 
of adequately large discrete systems. The general laws governing vortex-magnon interactions are 
established. We give simple physical explanations of the scattering results: the splitting of dou- 
blets for the modes with opposite signs of m, which takes place for the long wavelength limit, is an 
analogue of the Zeeman splitting in the effective magnetic field of the vortex. A singular behavior 
for the scattering amplitude, am oc k, takes place as k diverges; it corresponds to the generalized 
Levinson theorem and can be explained by the singular behavior of the effective magnetic field at 
the origin. 
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I. INTRODUCTION 

It is now firmly established that vortices play an im- 
portant role in the condensed matter physics of two- 
dimensional (2D) systems with continuously degenerate 
ground states. In particular, the presence of vortices in 
2D easy-plane (EP) magnets gives rise to the Berezin- 
skii-Kosterlitz-Thouless phase transitioniSi^ Vortices 
play an essential role in the thermal and dynamical prop- 
erties of 2D magnets, for a review see Ref. A vor- 
tex signature in dynamical response functions can be ob- 
served experimentally; e.g. translational motion of vor- 
tices leads to a central peak in the dynamical correlation 
functionsi^ This peak had been predicted predicted both 
by a vortex gas theory and by combined Monte Carlo - 
spin dynamics simulationsj^ 

Recently there has been renewed attention to the prob- 
lem of magnetic vortices for finite-size magnetic parti- 
cles, especially their dynamics. It becomes very impor- 
tant in connection with novel composite magnetic ma- 
terials such as magnetic dot arraysi&SiifliiLiSii^ These 
magnetic dots are submicron-sized islands made from 
soft magnetic materials on a nonmagnetic substrate. 
They are important from a practical standpoint as high- 
density magnetic storage devices, and are interesting 
as fundamentally new objects in the basic physics of 
magnetism. The distribution of magnetization in such 
a dot is quite nontrivial: when the dot size is above 
a critical value, an inhoniogeneous state with an out- 



of plane magnetic vortex occurs, which is stable due to 
competition between exchange and dipole interactionsii^ 
It is expected that these nonuniform states will drasti- 
cally change the dynamic and static properties of a dot 
in comparison with a uniformly magnetized magnetic 
disk. Recent experimentsi^iii verify such properties; in 
particular, a mode with anoma lously low frequency was 
detectedfii see also Refs. ll8ll9l 

The general properties of vortex dynamics are inti- 
mately connected to the problem of vortex-magnon in- 
teractions. Usually this problem has been studied numer- 
ically for discrete models, mainly for circular samples cut 
from large lattice systemsi^^i^-^i^^i^^i^'^i^^i^^i^''' An analyt- 
ical description of the problem in the framework of the 
continuum model has been proposed recently for differ- 
ent 2D magnets iS^iS2i22i2ai2i The most important effect 
of the vortex-magnon interaction is an excitation of cer- 
tain magnon modes due to vortex motion and vice versa. 
Because the magnons in the EP ferromagnet have a gap- 
less dispersion law, a possible Larmor dynamics of the 
vortex center is strongly coupled with a magnon cloud;^— 
therefore the corresponding motion has a non-Newtonian 
form.^^ 

In this paper we consider the magnon modes which 
exist in a 2D Heisenberg EP ferromagnet in the pres- 
ence of a vortex. We apply different, sometimes new, 
methods of analytical and numerical investigation, in 
order to extend the research of Ref. 0, presenting a 
wider range of results for the magnon scattering ampli- 
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tude. In Sec. m we demonstrate that the vortex acts on 
magnons in two ways. First of all, it provides for cou- 
pling between different directions of magnetization pre- 
cession; the magnon modes are described by a "general- 
ized" Schrodinger equation. Secondly, the problem nat- 
urally possesses an effective magnetic field, whose global 
properties are caused by the soliton topological charge. 
The scattering problem is treated both numerically and 
analytically for a wide range of wave numbers. The 
numerical study is carried out using two different ap- 
proaches: solving the eigenvalue problem for the contin- 
uum limit (in the weak EP anisotropy limit), and extract- 
ing the scattering data from numerical diagonalization of 
discrete systems, see Sec. lIIII The analytical study of the 
scattering problem is developed (Sec. llVjl using both the 
long- and short-wavelength approximations. In contrast 
to the previous worli24 we describe analytically the split- 
ting phenomenon of the doublets of magnon modes with 
opposite signs of the azimuthal quantum number, and 
give a physical picture of this effect: an effective magnetic 
field causes the Zeeman splitting of the magnon levels, see 
Sec. IIV Al The singular behavior of the scattering am- 
plitude is predicted in Sec. IIV Bl for the short wavelength 
limit. This feature is caused by the specific singular effec- 
tive magnetic field at the origin (vortex core); this study 
verifies the generalized version of the Levinson theorem, 
which we have established recently in Ref. for poten- 
tials with inverse square singularities. 



II. MODEL AND MAGNON MODES 

We consider the classical 2D Heisenberg ferromagnet 
(FM) with the Hamiltonian 



(1 - A)^^!^;,] , (1) 



where the spins Sn are classical vectors on a square lat- 
tice with the lattice constant a. Here (n, n') denotes 
nearest-neighbor lattice sites, J > is the exchange in- 
tegral, and A G [0, 1) describes easy-plane anisotropy. 

We consider the continuum dynamics of the Heisen- 
berg ferromagnet, which is adequate for model in the 
small-anisotropy case (1 — A^l). In a continuum limit 
the dynamics of the FM is described by Landau-Lifshitz 
equations for the normalized magnetization. 



m(r, t) 



S{r,t) 



sm o cos 0; sm t/ sm ( 



(2) 



The equations of motion result from the Lagrangian 



L = — / d\(l 



(3) 



with the energy functional 



E[0,, 



{Vef + {\7(j)fsin^e- 



(4) 



Here r^, = | •\/A/(l — A) is the characteristic length scale 
("magnetic length"). One should note that the strength 
of the weak easy-plane anisotropy (A ~ 1), determines 
the magnetic length scale r^; in this continuum analysis, 
results will depend on lengths scaled by r^, and have no 
explicit dependence on the anisotropy strength. 

The ground state of the magnet is continuously de- 
generate and isotropic within the easy plane (EP). The 
simplest elementary linear excitations of EP FMs that 
arise on the homogeneous background are the magnons 
belonging to a continuous spectrum. They have the form 
of cUiptically polarized waves; if one chooses the homoge- 
nous spin distribution along the x-axis, then the spin 
wave takes the form: 



= 1, mz+ imy 



The dispersion law has the gapless form 



oc exp(«A;a: — iut). (5) 



uj{k) = cky^l + k^rl 



(6) 



where c — 2aJS\/l — A is the characteristic magnon 
speed, k is the magnon wave vector, and k ~ \k\ is its 
magnitude. 

The simplest nonlinear excitation in the system is an 
out-of-plane (OP) vorte x??i?'^ 

<j) = <j)Q^ifo + qX, 6* = 0q{p), 

I-P (7) 
^o(0) = ^7r, eoi^) = TT/2, 

where ipQ is an arbitrary angle due to the EP symmetry, 
p = \r\/ry and x are dimensionless polar coordinates in 
the plane of the magnet, the vorticity g G Z plays the role 
of a TTi topological charge, and the polarization p = ±1 is 
connected with a tt2 topological charge (the Pontryagin 
index) : 



in 



Q — —Eij m ■ {dim x djm) . (8) 



We term Q the gyrocoupling density, following Thiele^; 
it has a sense as the density of the topological charge, see 
Papanicolaou and Tomaras'^^. For the vortex configura- 
tion ITJ the gyrocoupling density can be represented as 



Q 



q sm Oq ■ 



(9) 



hence the Pontryagin index takes on half-integer or in- 
teger values Q = qp/2. Note that the presence of a non- 
trivial 7r2-topological charge directly results in the gy- 
rotropical dynamics of the vortex, which conserves the 
gyrovector G — Q ■ 2TThSa~'^ez. 

The function 9o is the solution of an ordinary differen- 
tial equation, which can only be solved numericallyi^^iSi 
Without an external magnetic field two oppositely polar- 
ized vortices are energetically equivalent; for definiteness 
we set p = +1. 
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To analyze magnons on the vortex background, we use 
a formalism and set of coordinates developed in Ref. |^ 
setting up the problem in terms of local Cartesian spin 
components. The unperturbed spins of the static vor- 
tex structure, mo, define local polar axes 63, different at 
every site, specifically, So{r,t) — 6*63. 

It is to be understood that these axes depend on the 
site chosen. The magnetic fluctuations occur perpendic- 
ular to these local axes, suggesting the definition of other 
axes, 62 being chosen along the direction of x 63, and 
Ci = 62 X 63, to complete the mutually perpendicular 
set. One supposes that a dynamically fluctuating spin 
has small deviations along the ei and 62 axes so that a 
spin is written as 



S{r, t) ^ S (63 + miei + 771262) 



(10) 



The fields toi and 7712 have a simple physical significance, 
which can be seen if a given spin is supposed to have 
small deviations ip and i9 away from the vortex structure, 
determined by azimuthal and polar spherical angles, 0o 
and Oq. We write 



Sir,t) 
S 



= cos(0o + ^)ez + sin(6'o + i?) 
X [cos(0o + ^)ex + siii{(j)o + if>)ey 



(11) 



Linearizing in Lp and t?, and using the definitions of 
{ei, 62, 63}, comparison of Eqs. ifTUIl and lfTT|l shows that 



77ll = ??, 



7772 = LpSmOQ. 



(12) 



Thus, the mi field measures spin rotations moving to- 
wards the polar (e^) axis and the 7712 field measures spin 
rotations projected onto the xy-plane. In the absence of 
the vortex, we have 0o = 'I'/S, 4>o = 0, and such oscilla- 
tions correspond to the free magnons in the form jS)). 

The linearized equations for 7771 and 7?7.2 can be de- 
scribed by a single complex-valued function ip{r,t) = 
777,1 + *'772, which obeys the differential equation 

idr^ H^P + W^*, H={-iV-Af + U, (13) 

with the "potentials" 

Uip) = i sin^ 9o (1 - ^) - cos' ^0 - ^' (14a) 



mp)^isin^.o(l-^)H-f, 



A{p) 



q cos 6*0 



(14b) 
(14c) 



Here we use the dimensionless coordinate variable p — 
\r\/ry, dimensionless time variable t — t ■ c/ry, and the 
operator V = r^dr', prime denotes d/dp. 

Let us note that the vector A acts in the Schrodinger- 
like operator H in the same way as the vector-potential 
acts in the Hamiltonian of a charged particle. Then it is 



possible to conclude that there is an effective magnetic 
flux density 



B = V X A = e. 



q sin 9o ■ 6'q 



(15) 



Note that the effective magnetic flux density can eas- 
ily be rewritten through the gyrocoupling density J^J as 
B = Q ■ Ez- Therefore the total flux is determined by 
the nontrivial 7r2 topology of the vortex configuration. 
On first view when exploiting this analogy it is possi- 
ble to look for the Aharonov-Bohm phenomenon for the 
scattering problem, because this magnetic flux density is 
localized in the region of the vortex core. However, one 
can see that the total magnetic flux 



(16) 



is an integer multiple of the flux quantum $0 = 27r, so 
there is no Aharonov-Bohm scattering picture for the 
system. 

A differential equation like (|13|l is not a unique prop- 
erty of the vortex-magnon problem in the EP FM only. 
It appears for different kinds of anisotropy: it describes 
magnon modes on the soliton background in the easy- 
axis'^*' and isotropic magnetsSSi. Note that for the specific 
case of an isotropic system with an exact analytical soli- 
ton solution of the Belavin-Polyakov type, the potential 
W disappears, so the magnon modes satisfy the usual 
Schrodinger like equation idri' — H-tp, which describes, 
e.g., the quantum mechanical states for a charged particle 
in the axially symmetric potential U (p) under the action 
of an external magnetic field with a vector potential A. 

For the anisotropic case, when W 0, the prob- 
lem (|13|l has important unusual properties, which are 
absent for the Belavin-Polyakov case. More generally, 
there appear properties which are forbidden for the usual 
quantum mechanics. In particular, an effective discrete 
Hamiltonian of the system is not necessarily Hermitian; 
m Refs. l2lll2l25l some constructive methods were elabo- 
rated to avoid these problems. Nevertheless we will dis- 
cuss the features of Eq. H13|l in order to understand why 
the standard quantum mechanical intuition could fail. 

The standard quantum mechanical equation id-rip = 
Hip allows the conservation law StIV'P = ^ 3 for the 
current 



j = i {ipVip* - tp*Vip) + 2\iP\^A. 



(17) 



The generalized Schrodinger-like equation (|13l) with 
W violates this conservation law, namely: 



dr\lp\' 



-V j- iW 



(18) 



Nonconservation of probability density has posed some 
problems in the passage from standard quantum mechan- 
ics to old pre-Feynmann quantum electrodynamics. The 
reason is that Eq. H13|l is formulated neither for a Hermi- 
tian, nor a linear operator; the last statement is due to 
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the broken symmetry under the rescahng with 
A G C. There exists an analogy with relativistic theory: 
there can appear solutions with positive and negative en- 
ergy in the passage from the Klein-Gordon to Dirac equa- 
tion. In fact, our problem has the same origin. Let us 
reformulate the problem as an equation second order 
in time. One can calculate that the Klein-Gordon-like 
equation 



(19) 



is valid far from the vortex center. What is important is 
that Eq. (|19|) contains a Hermitian operator (similar ar- 
guments were used in Ref. I25I) . Therefore, the eigenvalue 
problem (EVP) for w^, not for uj, is more appropriate 
for this system; the only problem is to separate solutions 
with positive and negative u. Note, that there appear 
4*'^ order operators with respect to the space coordinate, 
which causes the presence of master and slave functions 
in the solution, see below. 

Such a problem, as well as a problem with non- 
conserved number of particles (probability amplitude), 
appears in the theory of a weakly nonideal Bose 
gas. It results, in fact, in the separation of positive 
and negative energy solutions under u-v Bogolyubov 
transformations . '^^ 

Following this scheme we need to generalize the u-v 
transformation to the nonhomogeneous case. 

We apply the partial-wave expansion, using the 
ansata^ 



V'(r,i)=E[^"(^)^''*°+«"('^)^"'*°]' 

a. 

$a(Xi t) =mx- UJat + = TUX - ^aT + IJ.n, 



(20) 



where a = (fc, m) is a full set of eigenvalues, m G Z 
being azimuthal quantum numbers, the rjm are arbitrary 
phases, and = ior^jc are dimensionless frequencies. 
This expansion leads to the following EVP for the radial 
eigenf unctions u and v (the index a will be omitted in 
the following): 



WW) = 12 H = 



W 



-W 



(21) 



+ 1X0 + 1/2 ± F is the 2D radial 



Here H± = 

Schrodinger-like operator with the "potentials 



U,{p)^U{p)+A^ + ^-\ 

2 

2m (A • e^) 2qmcos9o 



Vip) 



P 



\. (22) 
(23) 



Vp — jAp^ + (l/p)d/dp is the radial Laplace operator. 
In spite of the fact that the EVP H21|) is formulated for 
the Schrodinger operators iif±, this EVP is different in 



principle from the usual set of coupled Schrodinger equa- 
tions, which is widely used, e.g. for the description of 
multichannel scatteringi^ The reason is that the matrix 
Hamiltonian H is not Hermitian for the standard met- 
ric, for details see Ref. 25. To avoid this problem we 
introduce a corresponding bra-vector by the definition 



(24) 



The Hilbert space for the ?P'-function has an indefinite 
metric. 



= {u\u) - {v\v) 



(25) 



where {u\v) — u{p)v{p)pdp is the standard scalar 
product. By introducing such a Hermitian product, it 
is possible to define the standard energy functional, cp. 
Ref. US 

e[u,v] ^ {^\W\^) = iu\H+\u) + 2{u\W\v) + {v\H^\v). 

(26) 

Let us mention that Eq. (|21|l is invariant under the 
conjugations f2 —i- —Q, m — )■ —to, and u v. In a clas- 
sical theory we can choose either sign of the frequency; 
but in order to make contact with quantum theory, with 
a positive frequency and energy £fc = fiwfe, we will discuss 
the case !7 > (w > 0) only. Thus there are two different 
equations for the opposite signs of to. However, in the 
limiting case of the "zero modes" with J? = 0, the system 
again is invariant under conjugations to — > —to. For ex- 
ample, one of the zero modes, the so-called translational 
mode with to = +1, has the form 



sin Oq 



sin 0Q 
P 



+ ^0, (27) 



which describes the position shift of the soliton. Because 
of the degeneration of the EVP at i7 = 0, it leads to the 
existence of a zero mode with m = — 1; the eigenfunction 
of this mode can be expressed just from H27() under the 
conjugation u v. W e use here notations for mode 
indices as in Hefs. l29l8lt note that the mode with to = +1 
corresponds in our notation to the mode with to = — 1 
in the notations of Refs. 24 30, 

It should be stressed that the picture is quite different 
for the special case = 0, which corresponds to the 
isotropic magnet Here we have two uncoupled equa- 
tions for the functions u and v. One of the equations 
(for the eigenfunction v) has the negative eigenvalue —fi^ 
from which it necessarily results that w = 0. In this spe- 
cial case the zero modes (|27|l have the form u+i cx 9'^ and 
v^i — 0. Therefore the zero mode with to = — 1 cannot 
be obtained by the simple conjugation. It explains the 
difference between the collective dynamics of the soliton 
in isotropic magnets, where it is enough to take into ac- 
count only the mode with to = +1, and the EP FM, 
where translational modes with to = — 1 and m = +1 
must be taken into account. Nevertheless, the roles of 
the modes with to = — 1 and to = -1-1 are not equal, for 
details see Sec. IIV CI 
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III. SCATTERING PROBLEM: NUMERICAL 
RESULTS 

A. Continuum approach 

We intend to describe the scattering of niagnons by 
a vortex. However the EVP (|21|l is not adjusted for 
the scattering problem, because it does not provide the 
asymptotic independence of the equations at infinity. To 
solve the problem it is convenient to make a unitary 
transformation of the eigenvector | iP') , 



|a^)=A|4>), A 



cos £ — sm £ 
sin £ cos e 



(28) 



The angle £ of this unitary transformation is defined by 
the expression 



tan 2e = — . 

2Q 



(29) 



Then we obtain the following partial differential equation 
for the function 



H|^) = A|<^), H = Ho + V, A = di&g{K^-x^) , (30a) 
Ho = diag(J{o;-J^o), J{o--V2+lto, (30b) 
\/^[V + g-{W-l/2)]A-\ (30c) 



where g ~ || -i o || ^ metric spinor, the dimcnsionless 
wave number is k = fcr „ , and k = + 1. 

First let us consider the magnon spectrum in the ab- 
sence of a vortex (free fields). Without a vortex {q — 0, 

= 7r/2), Eqs. H30|l are uncoupled, which results in free 
magnons. 



Vm{p) = 0, 



2 / \m\TT TT 

• COS KO 7 

TTKp V 2 4 



(31) 



where Jm are Bessel functions. The free modes Um play 
the role of the partial cylinder waves of a plane spin wave 



exp {ik ■ r ~ iust) — i"^Jm{Kp)e 



imx—ii^t 



(32) 



with asymptotically independent solutions: 
C C 



Vp Vp 

1 ( ItoItt 

oc — — cos up — 

\fP \ 2 



TT 
4 



Vm\P) p 



9± 



->tp 



(34a) 
(34b) 

^m('«); 



The scattering results in the quantity 5„ 
it can be interpreted as the scattering phase shift, de- 
termining the intensity of the magnon scattering due to 
the presence of the vortex. Sometimes it is useful to 
introduce the scattering amplitude, Um = — tan (5m- Us- 
ing this notation, the oscillatory solution H34af) can be 
rewritten in the following form 



Um{p) oc J\m\{l^P) + CT„ • Y\„,\{Kp). 



3) 



where Y\ m\ are Neumann functions. Let us stress that the 
solution 134al|l is valid only in the sense of the asymptotic 
form H34a|) . As it follows from Eq. (|34bl) . the function v 
has an exponential behavior. 



Vm{p) OC K\^\{>tp) + 7„ 



{>tp) OC —— + 7, 



where K]^^\ and I\ 



343) 

are MacDonald and modified Bessel 
functions, respectively; at the same time u yields oscil- 
latory solutions. Naturally, the real modes have an os- 
cillatory form here; we will use this fact below for the 
numerical analysis. It means that the function u be- 
comes a master function in the Eq. H30a|) . while u is a 
slave (note, that we choose J7 > 0). This mirrors the dif- 
ference between Eq. (|30all and a usual set of Schrodinger 
equations. 

The scattering amplitude, or, equivalently, the phase 
shift, contains all information about the scattering pro- 
cesses. In particular, the general solution of the scat- 
tering problem for a plane wave can be expressed in the 
form, cp. Eq. |(3Jl 



m2 — imi- 



kVy 



oc e 



ik-r—iujt 



^iKp—iut 



where the scattering function has the form^S 

exp(-i7r/4) ^ , 2i5,„ . 



771— — OO 



m— — oo 



(35a) 



(35b) 



To describe magnon solutions in the presence of a vor- 
tex, one should note that far from the vortex center the 
potential V tends to zero, so the Eqs. H30a|l become un- 
coupled. 



(V^ + K^)u^ 0, (V^ - x^) u = 0, p > max ( 1; -; - 



(33) 



The total scattering cross section is given by the ex- 
pression 



E 



where q„ 
sections. 



(4/fc) sin^ (5m are the partial scattering cross 
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Let us switch to the numerical solution of the scatter- 
ing problem in the continuum approach. The differential 
problem to be integrated consists of Eq. H3()|l and asymp- 
totic conditions at the center of the vortex and at infinity: 



H|<i>) = A|«f), 



\^) - A 
1^) 



when p ^ 1 



when p 3> max ( 1 ; — ; — 



(36a) 
(36b) 

(36c) 



The presence of the matrix A in the condition Ij36b|) 
means that the functions u and v are not asymptotically 
independent even in the lowest approximation. In the 
next approximation there appears an additional "inter- 
action" between ii and which is realized in the nonunit 
factor em] its value cannot be found through this asymp- 
totic expansion. 

We use the one-parameter shooting method, solving 
Eqs. ijSni), as described in Ref. l3Qi>4Ck Choosing the 
shooting parameter €„, we "kill" the growing exponent 
for the function Vm m (j34bl|) . where the coefficient 7m 
should be equal to zero; as a result we have obtained a 
well-pronounced exponential decay for Vm oc K\m\i>^p), 
and oscillating solutions for u„i- The scattering ampli- 
tude was found from these data by comparison with the 
asymptotes (|36c|l . The results are discussed in Sec. lIII d 



B. Discrete approach 

In the discrete lattice approach, the small amplitude 
spin fluctuation modes in the presence of a vortex at the 
center of a finite circular system of radius R are found. 
The spins occupy sites on a square lattice. We use the for- 
malism and set of local coordinates as described in Sec.lITI 
for the continuum model. Similar to the continuum ex- 
pression (|10l) . we describe the dynamically fluctuating 
spin on lattice site n as 



Sn = S (63 + miBi + 771262) . 



cni) 



where mi and TO2 measure spin rotations moving towards 
the polar axis and projected onto the xy-plane, respec- 
tively, see Eq. ((T^ . 

The spin dynamics equations of motion with an as- 
sumed e~*"* time dependence were linearized in mi and 
TO2, leading to an eigenvalue problem requiring numer- 
ical diagonalization. We assumed a Dirichlet boundary 
condition, toi — TO2 = at the edge of the system stud- 
ied. For circular systems of radius R, we used a Gauss- 
Seidel relaxation scheme^^ to calculate the frequencies 
and eigenfunctions of some of the lowest eigenmodes with 
a single vortex present at the system center. Before do- 
ing this, the vortex structure was relaxed to an accurate 



static structure using an energy minimization scheme. 
The diagonalization is partial; typically only the lowest 
20 to 40 eigenstates were found, which substantially re- 
duces the computing time needed, and relaxes constraints 
on the precision of the calculations. This limited diag- 
onalization, however, gives only modes which have long 
wavelength spatial variations, which provides for a good 
comparison with continuum theory. 

We considered different values of A close to 1. Although 
the continuum limit would be better represented by us- 
ing A ver y close to 1, this could result in a vortex radius 

= I ^/X/{^^^^ easily exceeding the system size that 
can be treated numerically. Therefore, data were calcu- 
lated using A = 0.99, for which « 4.97a. With this 
size of vortex length scale, discreteness effects due to the 
underlying lattice should be unimportant, and still, the 
vortex structure fits well within the confines of a system 
with a radius as small as R w 10a, so that finite size 
effects should also be negligible. 

In general, a given mode has e™-^ angular dependence 
on the azimuthal coordinate Xi where m is some integer 
aziniuthal quantum number. In the continuum theory 
presented in Sec. m to is a good quantum number, due 
to rotational invariance. This symmetry is weakly broken 
on a lattice, but for long-wavelength and lower frequency 
modes, m can be considered a good quantum number 
even on a lattice. (Generally, the calculated magnon 
wavefunctions were sometimes found to be composed of 
linear combinations of +m and — m components.) The 
numerically found modes were also characterized by a 
principle quantum number n, being the number of nodes 
in the wavefunction along the radial direction. For a 
mode of determined m and n, the scattering amplitude 
a was found by a fitting procedure applied to the cal- 
culated eigenfunction for that mode (essentially, finding 
the ratio of outgoing and incoming waves), see Ref. 1251 
for details. 

In the continuum theory, scattering was analyzed as a 
function of wavevector k, or in terms of the dimension- 
less kvy. For lattice calculations, the values of k cannot 
be chosen freely, instead, they are determined by the ac- 
tual system size. For a mode found to be oscillating at 
eigenfrequency w, the wavevector magnitude k associated 
with the mode was found by supposing k = (fc, 0), and 
inverting the free magnon dispersion relation for the 2D 
EP FM on a lattice, 



4J5V(l-7fe)(l- A7fc), 



7fc = -^{cosk^ + cos ky). 



(37) 



Therefore, a calculation of the modes for a single lat- 
tice size gives only specific values of k = kr^, one value 
corresponding to each mode. To get a wider and more 
continuous range of data for comparison with the con- 
tinuum theory, calculations were carried out on lattices 
ranging in radius from R ~ 15a to i? = 40a. By plotting 
results as functions of kr^ , for fixed m but from various n 
and _R, the data from the different system sizes superim- 
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FIG. 1: (Color online) Scattering data for different m 
for small wave numbers, kvv < 1.3: from continuum the- 
ory (lines) and from discrete model numerical diagonaliza- 
tion (symbols) in circular square lattice systems of radii 
R = 15,20, 25,30,35,40. 



poses smoothly, giving more slowly changing fcr^,, which 
is more appropriate for comparison with the continuum 
limit. 



FIG. 2: (Color online) Scattering data for different m for a 
wide region of wave numbers k: from continuum theory (lines) 
and from discrete model numerical diagonalization (symbols). 



there is no real divergence; it means that the physically 
observed phase shift does not tend to zero at infinity, 
but to a finite value Sm{k 00) ±7r/2. The scatter- 
ing data are presented in Figs.n]|5| Comparison with the 
results of exact diagonalization on finite systems shows 
very good agreement between the two approaches. 



Numerical results 



Numerically, we have obtained the data of the vortex- 
magnon scattering by the two different approaches dis- 
cussed above: solving the scattering problem using 
the shooting method for the continuum limit, and ex- 
tracting the scattering data from numerical diagonaliza- 
tion of finite discrete systems. To be specific, data are 
presented for scattering from a vortex with unit vorticity, 
^ = -1-1 and positive polarization, p = +1. One should 
note that results for vortex-magnon scattering for modes 
m from other vortex types, as seen in Eq. H23|l . should be 
depend on the sign of qpm. The results arc the following: 

For all modes the scattering amplitude am{k) tends to 
zero as — > 0. In the long- wavelength limit the maximal 
scattering is related to the modes with m = ±1. Except 
for the mode with m = — 1, the scattering amplitude 
in the long-wavelength limit takes a negative value, see 
Fig.n At extremely low values of wave number k < 0.01, 
the scattering data contain sets of doublets for modes 
with opposite signs of m. In the long-wavelength limit 
the doublet splitting appears as a small correction, but 
the scattering picture changes when k increases. For all 
modes am{k) diverges as fc — > 00: the scattering ampli- 
tude am{k 00) — !■ -l-oo for all modes with m > —1, but 
Cmik — > 00) —00 for TO < —1, see Fig. |21 Naturally, 



IV. SCATTERING PROBLEM: ANALYTICAL 
DESCRIPTION 

A. Scattering at long wavelength 

In order to analyze the scattering problem analytically 
in the long-wavelength limit, we start from the zero- 
frequency solutions, when i7 = 0. First note that for 
the special cases to = 0,±1 there exist so-called half- 
bound states. Recall that a zero-frequency solution of the 
Schrodinger-like equation is called a half bound state if 
its wave function is finite, but does not decay fast enough 
at infinity to be square integrable. We will refer to such 
modes as half-local modes. These modes correspond to 
the translational (m — ±1) and rotational (to = 0) sym- 
metry of an infinite system, they have an exact analytical 
form: 

-^^ = ^^'^^"^=-■^0, m = 0,±l. (38) 

Unlike the case of half-local modes with m = 0,±I, all 
other zero-frequency solutions are nonlocal, and we are 
not able to construct exact expressions for them, but only 
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the asymptotes for p 1: 



the ansatz 



(39) 



Nevertheless, we will see that the knowledge of asymp- 
totic solutions like (|39|l will be enough to reconstruct 
the K-dependence of the scattering amplitude. In order 
to solve the scattering problem in the long wavelength 
limit we apply a special perturbation theory, proposed 
in Ref. for the modes with m = ±1,0, and extend- 

ing it for all values of to. We construct the asymptotes of 
such a solution for a small but finite frequency by making 



v{p) = i'o(p) ■ H- Kai(p) + k^q;2(p) , (40a) 



u{p) = Mo(p) ■ 1 + K/3i(p) + K^/32(p) 



(40b) 



here ol\, ji\ and a2, /32 are first and second order correc- 
tions to the zeroth-solutions, respectively. Let us insert 
this ansatz into the set of Eqs. H36|) . multiply from the 
left with p • without integrating; then one obtains 
equations for the first and second order corrections: 



1,2, 



$2(p) = p\ ~ul + vll2 + 2{W - 1/2) {ul - «2) + 2V {ulP^ - vla^) + lu^v^ {{W - 1/2) (ai + /3i) - 2V 



(41) 



We are interested in the corrections which will give us a possibility to calculate the scattering amplitude. The 
formal solution of these equations can be written as 



/3fc(p)=/3(0) 



P afc( ?7)«o(?/)d?7 
ulij]) 



dr] 



wliv) Jo 



Mm- 



(42) 



r 



Let us calculate the first order correction (3i . It is easy to 
see that the second RHS-term has an exponential decay 
as p ^ oo, while the third one has a slow algebraic decay 
only. Thus, far from the vortex core we have simply 



/?! (p) ~ const 



0-2|m| 

2|to| 



(43) 



valid in the region p 1. 

To calculate the second order correction, /?2, let us 
note that the last RHS-term of Eq. H42|l is divergent for 
p oo, while the integral with a2 has an exponential 
decay, like the first order correction. To derive the diver- 
gent inner integral in Eq. (|42|l we add and subtract the 
function 



[sin^ eoul]' 



2(|to| + 1) 



2 TO 



(44) 



Then we arrive at an approximation for /32 (p) in the im- 
portant region p ^ 1: 



P2{p) — const 



-2|m| 



4(|r 



2 TO 



Inp 
2M 



(45) 



Now we are in position to compare the magnon am- 
plitude Um = uo{l + K/3i -I- K^/32) with the scattering 



approach in order to extract the information about the 
scattering amplitude amin). To describe the scatter- 
ing problem in the long-wavelength approximation we 
rewrite the differential problem H30|) for large distances 
Kp ^ 1, only considering the terms with k^. In this scat- 
tering approach the oscillating function iim satisfies an 
equation 



0, 



The solution of this equation can be written as 

(p) cx J\^\{kp) + a,y{K)Y\^i{Kp) 
1 f W\n 



cx COS 



np - 



2 4+'" 



(46) 



(47) 



where the index of the Bessel and the Neumann function 
is noninteger. It results in a value of dm which differs 
from the real scattering phase shift Sm- Using asymptotic 
expansions (|34a|) and 147() . the desired relation between 
the phase shift and Sm can be written as 



(5m (k) = S^{k) + 



(48) 



In the lowest order approximation in k, the corresponding 
relation for the scattering amplitudes has the form: 



4 TO 



(49) 
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To compare the scattering solution (|47|l with the result 
of the perturbation theory we can expand the cylindri- 
cal functions in powers of the small quantity \i>\ — \m\ 
and represent through the cylindrical functions of inte- 
ger order |m|, as done in Ref. 31. After that in the region 
Kp ^ 1 we are able to use the asymptotes of the cylin- 
drical functions at the origin; we arrive at the formula 



- 1 1 - , - In f ^) 

^ I 4(|m| + 1) 2\m\ \ 2 J 

(|to|!)2 / 2 yi"! r 



7r|m| \Kp 

|m|-l 



1 - 



4|m| 



>, (50) 



Srn =J+\m\ ^ 



n{\m\ — n) ' 



where 7 is Euler's constant. 

Comparing this expression with the perturbation the- 
ory results [see Eqs. (|40b|l . (|^ ]. in the region 
1 p ^ 1/k, where both are valid, we can restore the 
general dependence of the scattering amplitude in the 
long-wavelength approximation: 



. 2-K\m\ 



2\m\ /^\2|m| + l 

mT>rn. ( - I , (51a 



5n(|m|!)2 J, 



m(|m|!)2 7o 



<i'2(e)-'i>f (o]de, (51b) 

$i(Ode (51c) 



Eq. (|51a|l solves the scattering problem except for factors 
Am and S^. These factors can be found by the numeri- 
cal integration of Eqs. (jSlbp and (|51c|) . using numerical 
data for uq and vq. Thus, solving the equations for zero- 
modes once, we compute the whole dependence crm('t). 
Nevertheless, in order to discuss the analytical behavior 
let us note that for sufRciently large values of |m|, we can 
limit ourselves to the contribution of the term with Mq in 
the function $1 and the term with Uq/3i in the function 
$2, see Eq. (|^ . To calculate the integrals we need to 
have more information about the zero-modes. At small 
distances p <C 1 the isotropic (exchange) approximation 
works correctly, which leads to the following solutionsi2& 



uq oc p'™' sin 6*0, vq cc p 



\m\ + lQ, 



Such solutions have the correct asymptotic behavior at 
infinity and at the origin. Our numerical calculations 
justify the correctness of these assumptions for m ^ 1; as 
a result we obtain analytical estimates for these factors: 



18(2|m| - 1)! 



'Br. 



47r(2|TO| - 1)! 



(HO 



n2 



(52) 



taking into account 151a|) . and asymptotes for half-local 
modes from Refs. I24I41I 



o-m=o(K) « ^2'^^ ln(l/K), 

TTK 

0-m=±l(K) ~ T — , 



(53a) 
(53b) 



/ f^\2\m\ /^\2|m| + l 

O'm#0,±l(K) ~ --^m (^2 j +m'B,n(^-j ■ (53c) 



In the main approximation in k the scattering picture 
contains doublets for modes with opposite signs of m 
for the modes with \m\ > 1. The splitting of the dou- 
blets (the last term in Eq. (|53c|) 'l appears in the next 
order in k. The splitting of the doublets for the magnon 
modes on the vortex background with given m = ±n 
was mentioned in the earliest papers on vortex-magnon 
scattering)2iiS^ but it was explained, in fact, only for 
m = ±li24 Our considerations on the basis of Eq. H13() 
show that the splitting of the scattering data is the direct 
analogue of the Zeeman effect for electron states splitting 
in an external magnetic field. To follow this analogy one 
can rewrite the splitting constant in the form: 



de e'l'"! {AiO ■ e. 



hence the splitting appears only in the effective magnetic 
field, which is described by the vector potential A. 

Using scattering results H53|l one can solve the scatter- 
ing problem for a plane spin wave in the form (|35|l . In 
the long-wavelength limit the maximum scattering is re- 
lated to the translation modes with m = ±1, which gives 
the scattering function l|35bp in the form 



e^^/^inx. 



(54) 



In this approximation the scattering is anisotropic, and 
the total scattering cross section is g = tt'^k/A. To explain 
the origin of the anisotropic scattering, let us mention 
that the plane spin wave makes a spin flux, which influ- 
ences the vortex as a whole, trying to move it by exciting 
translational modes. It is well-known that the vortex dy- 
namics appears in the gradient of a magnetization field 
like the magnon fluxiiSi The dynamics of the vortex has a 
gyroscopical behavior (see for the reviewQ): acting along 
the X-axis, the spin wave causes the translational motion 
of the vortex along the y-axis, which results in H54|l . 



B. Scattering problem for short wavelength 



To compare the scattering results for different modes, we 
write explicitly the asymptotic expressions for all modes. 



For large fc, in the main approximation to lowest order 
in 1/k, the scattering problem H3()|) can be rewritten in 
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FIG. 3: (Color online) Scattering data for different m for the 
short -wavelength limit: from asymptotes l|67|l of the contin- 
uum theory (lines) and from discrete model numerical diago- 
nalization (symbols) 



the form: 


H^) = 











„|?n — 1 



when p <§; 1 



K\m\{np) 



(55a) 
(55b) 

when p ^ 1. 

(55c) 



We see that the functions u and v have independent 
asymptotes at the origin (|55b|l and at infinity H55c|l . It 
means that the role of the "coupling potential" W in the 
scattering problem (|55|l is unimportant here. Therefore 
one can neglect the "coupling potential" and formulate 
the scattering problem for the master function u only: 



where the partial potential is 



(56) 



U^{p)^Uo{p) + V{p)^U{p)- 



m 
P 



(57) 

It is natural to suppose that the WKB-approximation 
is valid for this case. We use the WKB-method in the 
form proposed earlier for the description of the scattering 
for isotropic 2D magnetsSS, and generalized after that 
for any singular potentials. ^'^ We start from the effective 
ID Schrodinger equation for the radial function Um{p) = 



'4'm{p)/^/p, which yields 



— + 



dp- 

Woff(p) = %n{p) 



1 



(58) 



The WKB-solution of Eq. f^ . i.e. the ID wave function 
V'm^'^j leads to the following form of the partial wave 



-WKB 



WKB 



VP 



-r=r^ cos f xo + np')dp') , 

(59) 



where T = ^/K^~^l^s. Analysis shows that Eq. itS^ is 
valid for p > a, where a is the turning point. The value 
of a corresponds to the condition CP(a) = 0, which results 
in a ~ \m\/K <^ 1. We assume that the parameter po 
satisfies the condition a <ti pa ^ 1- 

On the other hand, at small distances p ^ 1, the par- 
tial potential Um has the asymptotic form 



If ^ 

P^ 

V — m — lim \p {A{p) ■ e,,)] = m + qp, 



(60) 



therefore one can construct asymptotically exact solu- 
tions (recall that we suppose q = p = 1) 

Um oc J|„+i|(kp), when p < 1. (61) 

For K ^ |rn| there is a wide range of values of p, namely 

|m|/A:<p<l, (62) 

where we can use the asymptotic expression?^, for the 
Bessel function (|61|l in the limit kp ^ |to|: 



1 / to + Itt TT 4m+l2-l\ 

Umoc— cos Kp -H . (63) 

^ V 2 4 



8kP 



In the range of Eq. H62|l the solutions H59|) and (|63|l coin- 
cide due to the overlap of the entire range of parameters, 
so one can derive the phase xo in the WKB-solution H59|l. 



\m + 1 1 TT 



Xo = Kpo 



IT 

4 



4|to-|- 1|2 - 1 
8kpo 



Therefore, we are able to calculate the short-wavelength 
asymptotic expression for the scattered wave phase shift 
by the asymptotic expansion of the WKB-solution l|59|l : 




+ 



■?{p')dp' + Xo - up 



ImlTT TT 4to^ 



(64) 



8kp 
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Under assumed conditions kp^ 1, the WKB-integral in 
(|(j4|l can be calculated in the leading approximation in 
1/kp, 



y{p)dp' w k{p - Pa) 



2k 



Ucs{p )dp . 



As a result, the scattering phase shift for large wave num- 
bers, \, has the form 



(5„(k) = 5,n{oo) ~ I Um{p) - ^ 



2k 



dp, 



with the limiting value 

l5m(oo) : 

sgn_|_(TO) = 




--•sgn+(m), (65) 



Calculating the integral we obtain the phase shift in the 
form: 



(>m[K) = ■ sgn+(m) + 



(66) 



1 r f 3 sin^ ^ g ^^^2 ^ ^2 ^ ^ 2.44, 



2), 



k 

1 — cos 6*0 



P 



dp w 1.38. 



The corresponding amplitude of the vortex-magnon scat- 
tering is 



Hi -f mT>2 



(67) 



This linear divergence is well-pronounced in the numer- 
ical results, see Fig. |31 To understand the origin of this 
divergence let us go back to the Eq. I|65|) . One can see 
that the scattering phase shift at fc ^ cxd does not van- 
ish for potentials with inverse square singularity at the 
origin, with v ^ see Eq. ()6U|I . This is possible only 
in the magnetic field, which has singular behavior like 
\A\ ~ l/p. 

Let us look for the consequences of this unusual behav- 
ior of the scattering, a„i ±oo. We consider the scat- 
tering problem for a plane spin wave in the form (|35|1 . In 
the short-wavelength limit the WKB results for the phase 
shift (|66|l are available. One can see that the scattering 
function (j35b|l tends to zero very quickly for large wave 
numbers, 9^(x) = 0(k~^/^), so there is no real divergence 
or singularity for a physically observable quantity such 
as the total scattering function "J at large energies. 



C. Levinson theorem 

Now we can compare the scattering results in the long- 
and short-wavelength limits. The scattering is absent 



is 



CD 



CO 

o 
C/5 



'^^ -0.1 



CD 
(/> 
CO 



CO 
o 



-0.2 



-0.3 



-0.4 



-0.5 



(a) m<-1 




1 1.5 2 2.5 
Wave number, kr^ 

(b) m>-1 



m=-1 

m=0 

m=1 
m=2 

m=3 

IT1=4 

0.5 1 1.5 2 2.5 3 
Wave number, kr^ 



3.5 



FIG. 4: (Color online) Scattering phase shifts for different 
m. Numerical results from the continuum theory. 



for the limit fc — )■ 0. However, the scattering amplitude 
has a linear divergence a oc k for sufficiently large wave 
numbers, see Eq. (|67|) . All these results were verified 
by the numerical calculations for continuum limit and 
for finite sized discrete lattice systems, see Figs. ^ |3 
O According to our analytical calculations, see Eq. H65|l . 
the phase shift for the short-wavelength limit tends to 
the finite value Sm{oo) = — sgn^(m) • 7r/2. This result 
corresponds to the numerical data, see Fig.01 except for 
the mode with m = — 1, where the numerical data gives 
Sm=-iioo) = —tt/2. However we need to note that the 
phase shift is determined with respect to tt, in this sense 
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values Sm = tt/2 and 5m = —tt/2 are identical. What 
is physically important is how Smi^) changes from small 
to large k. According to our numerical results, the total 
phase shift can be described by the formula 

<5m(0) - d„i(oo) = ^ ^ (68) 
[2, m=-l. 

It is well-known that the total phase shift is related to 
the number of bound states according to the Levin- 
son theorem from the scattering problem for a spinless 
quantum mechanical particle without a magnetic field. 
This theorem was originally proved by Levinson for the 
3D case, see the review nSft. The two-dimensional version 
of the Levinson theorem reads^^i^^i^ 

6rn{0) - dmicx^) ^ n ■ (69a) 

If there exist half-bound states (see notations in the 
Sec. lIV A^l for the p-wave (m = 1), this is modified to^^i^ 

6i (0) - (5i (oo) = TT • iV|= + TT. (69b) 

For the 2D EP FM, the scattering picture is much more 
complicated. First, we have no standard Schrodinger 
equation, but the generahzed one l(T^ . This becomes 
apparent at most in the threshold behavior for the half- 
bound states, and the contribution of the half-bound 
states in the form Ij69bp may be not adequate, see be- 
low. Second, because of the role of the effective magnetic 
field, there appears an m-dependent potential: the sym- 
metry Smin) = S-m{K) is broken, so it is not enough to 
take into account partial waves with m > only. As a 
result Levinson's relation (|69a|l has a different form for 
the opposite signs of m. 

Thus, except for the case of half-bound states one can 
hope that the Levinson theorem is adequate. However we 
see that the total phase shift H68|l contradicts the Levin- 
son theorem in the form H69I) . The reason is that the 
partial potential ILm in the Schrodinger equation (|56|l has 
an inverse square singularity at the origin, ILm ^ ^'^ I ^ 
where v — •m-\- qp, see ()60(l . Such a situation changes the 
statement of the Levinson theorem. As we have proved 
recently in Ref . ^2, the generalized Levinson theorem for 
the Schrodinger-like equation for potentials with such 
singularities has form: 

5m{Q) - Sm{^) = ^ • A^„U f (IH - H) ■ (70) 

An additional tt can appear on the RHS of this equation, 
if the half-bound states exist for the p-wave {\m\ — 1), 
see (j69bll . To explain the meaning of the extra term 
(7r/2) • (|j/| — |to|) in the generalized Levinson theorem 
(|70|l . recall that in the partial wave method the scatter- 
ing data are classified by the azimuthal quantum number 
TO, which is the strength of the centrifugal potential. In 
the presence of a partial potential with an inverse square 



singularity at the origin such as Um ~ I the effec- 
tive singularity strength is shifted by the value \v\ — |to|, 
which results in a change in the short-wavelength scat- 
tering phase shift by (7r/2) • (|to| — \v\)- 

Let us compare the predictions of the general- 
ized Levinson theorem (|70|l . which is suitable for the 
Schrodinger-like equation, with our results for the 
vortex-magnon scattering problem in the 2D EP FM, 
which can be described by the generalized Schrodinger 
equation p3|l. In our case the singular potential is 
caused by the specific singular magnetic field at the ori- 
gin, \A\ ^ 1/ p, which results in v = \m + 1|. The system 
has no bound states, Nm = 0, therefore Eq. (|70|l takes 
the form: 

TT 

Smifi) - 5m{.O0) = - ■ Sgn_|_ TO. (71) 

Our numerical results H68(l correspond to this formula for 
all modes with m 7^ — 1. The cause is the influence of the 
half-bound states. By comparison of H71|l and (|68|l . one 
can adapt the generalized Levinson theorem for this case. 
It reads 

5m{^) - Smico) 

^{n-N^ + ^{\v\-\m\), when -I (72) 

~ y-K ■ N^^ + TT +^ - \m\) , when to = -1. 

Let us compare this result with Eq. H69|) . An extra tt, 
which appears for the mode to = —1, is connected with 
the half-bound states, see Eq. Ij69b|l . To explain the 
situation, let us stress again that our scattering prob- 
lem is formulated not for the standard Schrodinger equa- 
tion. However, the problem has a symmetry such that 
one eigenfunction becomes a master function, while the 
other is a slave. This makes it possible to use the main 
features of the standard quantum mechanical scattering 
theory. The appearance of the half-bound states is con- 
nected with the symmetry of the whole system, and both 
of the eigenfunctions are important. In the system there 
are three half-local modes, see Eq. (|38|l . According to 
Eq. (|72|l only one of the half-bound modes, namely, the 
mode with to = —1, gives an extra tt to the Levin- 
son's relation. More generally, this extra contribution 
corresponds to the half-bound mode with to = —qp, 
see Eq. (jBUJ. This result cannot be explained in the 
framework of the Levinson theorem for the standard 
Schrodinger equation, where both half-bound states with 
TO = -|-1, and TO = — 1 should make contributions to 
the Levinson's relation. The corresponding analogue of 
the Levinson theorem for the generalized Schrodinger 
equation H13|) takes into account the contribution of the 
half-bound state for only one value of to, namely, for 
TO — ~qp. 

V. CONCLUSION 

We have presented a detailed study of vortex-magnon 
interactions in the 2D EP ferromagnet, having described 
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this process by a "generalized" Schrodinger equation. 
The main features of the magnon scattering are con- 
nected with the special role of the effective magnetic 
field, which is created by the vortex. This effective field 
acts on magnons in the same way as a magnetic field 
influences an electron, leading to the appearance of the 
Lorentz force and the Zeeman splitting of the magnon 
states with opposite values of the azimuthal numbers m. 
The singular behavior of the effective magnetic field at 
the origin causes a divergence of the scattering ampli- 
tudes for all the partial waves; we have confirmed this 
study by a generalized version of the Levinson theorem 
for potentials with inverse square singularities. 

Our investigations can be applied to the description 
of the internal dynamics of vortex state magnetic dots; 
the theory of the vortex-magnon scattering developed 
here could be a good guide for the study of the normal 
modes in vortex-state magnetic dots. It is clear that 
the EP FM cannot correspond quantitatively to the case 
of vortex-state magnetic dots, where the anisotropy is 
negligible and the static vortex structure is stabilized by 
magnetic-dipole interactions. We did not consider this 
type of interaction in this paper, as it is difficult to ac- 



count for. Nevertheless, we believe that the main features 
of the problem studied above are generic. For example, 
an effective magnetic field exists due to the topological 
properties of the vortex only. Furthermore, we expect the 
appearance of modes with anomalously small frequencies, 
e.g., the mode of the translational oscillations of the vor- 
tex center. The nonzero frequency of this mode is caused 
by the interaction with the boundary only. Additionally, 
the splitting of the doublets for modes with opposite m 
should appear due to the role of the effective magnetic 
field. 
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